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We propose a general theory for the critical and pseudogap temperatures Tc and T* dependence 
on the doping concentration for high- Tc oxides, taking into account the charge inhomogeneities in 
the Cu02 planes. Several recent experiments have revealed that the charge density p in a given 
compound (mostly underdoped) is intrinsic inhomogeneous with large spatial variations which leads 
to a local charge density p(r). These differences in the local charge concentration yield insulator and 
metallic regions, either in an intrinsic granular or in a stripe morphology. In the metallic region, 
the inhomogeneous charge density produces also spatial or local distributions which form Cooper 
pairs at a local superconducting critical temperatures Tc{r) and zero temperature gap Ao(r). For 
a given compound, the measured onset of vanishing gap temperature is identified as the pseudogap 
temperature, that is, T*, which is the maximum of all Tdr). Below T* , due to the distribution 
of rc(r)'s, there are some superconducting regions surrounded by insulator or metallic medium. 
The transition to a coherent superconducting state corresponds to the percolation threshold among 
the superconducting regions with different Tc(r)'s. The charge inhomogeneities have been studied 
by recent STM/S experiments which provided a model for our phenomenological distribution. To 
make definite calculations and compare with the experimental results, we derive phase diagrams 
for the BSCO, LSCO and YBCO families, with a mean field theory for superconductivity using an 
extended Hubbard Hamiltonian. We show also that this novel approach provides new insights on 
several experimental features of high- Tc oxides. 

Pacs Numbers:74.72.-h, 74.20.-z, 74.80.-g, 71.38.-fi 



I. INTRODUCTION 

It is well known that the properties of high- 
temperature superconductors (HTSC) vary in an unusual 
way when a moderate density of holes are introduced into 
the Cu02 planes by chemical doping. This is one of the 
reasons why, despite of a large experimental and theoret- 
ical effort, the nature of the superconductivity in these 
materials remains to be explained'^. Correlations make 
the parent undoped compound to be a Mott insulator 
and, upon doping, the underdoped compounds display 
unusual metallic properties with increasing Tc- Doping 
beyond the optimal level yields normal metals with Fermi 
liquid behavior and with decreasing Tc. 

The non usual properties of underdoped samples have 
motivated several experiments and two features have 
been discovered which distinguish them from the over- 
doped compounds: firstly, the appearance of a pseudogap 
at a temperature T*, that is, a discrete structure of the 
energy spectrum above Tc, identified by several different 
probes^. T* is found to be present also in overdoped 
samples^''^, but at temperatures near Tc. Second, there 
are increasing evidences that the electrical charges are 
highly inhomogeneous up to (and even further) the opti- 
mally doped region"*"^ . These charge inhomogeneities are 
neither due to impurities nor to crystal defects, but are 
intrinsic to the type of cuprate, producing local lattice 
distortions in the Cu02 bond length^". In fact, such in- 
trinsic inhomogeneities are also consistent with the pres- 
ence of charge domains either in a granular^"® or in a 



stripe^^ ^'^ form. 

Therefore, in our view, it is very likely that the pseudo- 
gap and the intrinsic charge inhomogeneities are closely 
related and understanding their interplay is of great im- 
portance to understand the general phenomenology and 
the phase diagrams of the HTSC families. Currently 
there are two main different proposals to explain the ex- 
istence of the pseudogap: In the first one, the pseudogap 
is regarded as a normal state precursor of the supercon- 
ducting gap due to local dynamic pairing correlations in a 
state without long range phase coherence, and Tc is much 
smaller than T* because of strong phase fluctuations^'*' 
In the second proposal, the pseudogap is a normal state 
gap, which is necessarily independent of the supercon- 
ducting gap and which competes with the superconduc- 
tivity, existing even below Tc for compounds around the 
optimum doping, ending in a quantum critical point^^^^^ 
at zero temperature. 

Recently we have described a simple new scenario^^: 
when the temperature of a given sample is decreasing 
and reaches T*, it activates the formation of pairs at 
some selected low doping metallic regions. Initially these 
superconducting or pair-rich regions are not connected 
and there is no phase coherence over the whole sys- 
tem. The coherent superconducting transition occurs 
when the temperature reaches a value (Tc) at which the 
different superconducting regions percolate. This sce- 
nario relies heavily upon the fact that the cuprates have 
an intrinsic inhomogeneous charge distribution. Thus 
a given HTSC compound with an average hole per Cu 
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ion density (p) and with an inhomogeneous microscopic 
charge distribution p{r) has a distribution of small clus- 
ters or stripes with a given local Tc{r) (0 < (p) < 1 
and < p{r) < 1). Thus the model depends strongly 
on the intrinsic charge distribution p{r) although their 
exactly form is not well known and since the discovery 
of the spin-charge stripes^^, they are a matter of intense 
current research. As a consequence, several recent exper- 
iments have demonstrated that the charge distributions 
are more inhomogeneous for underdoped and more ho- 
mogeneous for overdoped compounds and may resemble 
either a granular^"^ or a stripe structure. 

In the spin-charge striped scenario, some regions of the 
plane are heavily doped (the stripes) and other regions 
are underdoped, filling the space between the charge-rich 
stripes. Stripe phases occur due to the antiferromag- 
netic interaction among magnetic ions and Coulomb in- 
teraction between the charges, both of which favor lo- 
calization. On the other hand, the zero-point motion of 
the holes favors delocalization and tends to create phase- 
separated states rich either of spin or charge. Experimen- 
tally, the stripes are more easily detected in insulating 
materials, where they are static, but there are evidences 
of fluctuating stripe correlations in metallic and super- 
conducting compounds^^'^^'^'^. 

Based on the results of the above mentioned exper- 
iments, we have introduced a local charge distribution 
p{r) to model the real charge distributions inside a HTSC 
compound. Since doping produces an intrinsic inhomo- 
geneous charge distribution and metallic and insulator 
regions seems to coexist, we assume that it contains two 
parts: one for a hole-rich and other for a hole-poor par- 
tition. This mimics the striped phases. The hole-poor 
regions are, in most cases, antiferromagnetic Mott insula- 
tors. The hole-rich regions form an inhomogeneous metal 
with spatially varying charge density. The existence of 
the intrinsic charge inhomogeneities have several conse- 
quences, and one of them is the non-Fermi liquid behavior 
of the underdoped compounds. As concerns the super- 
conductivity, they produce spatially dependent supercon- 
ducting gaps Asc{r) due to short coherence length and 
also spatially dependent superconducting domains with 
critical temperatures Tc{r). Recently, very fine scanning 
tunneling microscopy/spectroscopy (STM/S) data^'^ has 
revealed the spatial variation through the differential 
conductance which provides a strong evidence for such 
distribution of zero temperature superconducting gap 
/:^sc{r)- More recently, new STM data on Bi2212 us- 
ing scattering resonances at Ni impurities atoms, has re- 
vealed a large nanoscale spatial dependence of the super- 
conducting gap^. 

All together, these experiments show variations on lo- 
cal density of states (LDOS) of just a few angstroms. 
This is consistent with the very short coherence length in 
HTSC. Since T'c(r) is proportional to Asc(?'), the open- 
ing of the largest gap occurs at the highest of all the 
Tc(r)'s, which is exactly T* . Above T* there is no gap 
and below it, there is the development of some super- 



conducting clusters inside the material. Depending on 
the distribution of the values of Tc(r) in the compound, 
some regions become superconducting. Decreasing the 
temperature, the number of superconducting clusters in- 
creases, bringing about the superconducting domains to 
grow. When the temperature reaches a value which is the 
Tc of the compound, the superconducting regions per- 
colate through the sample and therefore it can hold a 
dissipationless current. 

On Fig.l we show a schematic local charge density dis- 
tribution p{r) which is based on the above experimental 
data, as mentioned, and could give rise to a gap dis- 
tribution of Asc{f)- This is a pictorial sketch of what 
should be the real charge distribution for a optimal doped 
HTSC, based on the above information. It reflects the 
nano-scale variations of the superconducting energy gap 
or the nano-scale variations on the density of charges in 
the material. 




FIG. 1. Sketch of the possible charge distribution in a 
HTSC compound. The white part represents the insulator 
or hole-poor regions. The colored region mimics the hole-rich 
metallic part of the compound. The lighter colors are for 
metals with higher values of the superconducting gap and the 
darker regions is for lower gap regions. 

This percolating scenario can be understood by analyz- 
ing the scanning SQUID microscopy magnetic data which 
makes a map of the expelled magnetic flux (Meissner ef- 
fect) domains on LSCO films^". This experiment shows 
the regions where the Meissner effect continuously devel- 
ops from near T* to temperatures well below the perco- 
lating threshold Tc- Notice that the regions without mag- 
netic flux due to the Meissner effect can exist only where 
the Cooper pairs are in phase coherence in a large spa- 
tial region. Furthermore there is no apparent difference, 
other than the size of the superconducting region, as the 
temperature changes across Tc. Secondly, it shows the 
existence of the superconducting regions above Tc which 
differ only in size from those below Tc. Furthermore, the 
d-wave nature of both pseudogap and order parameters, 
the tunneling conductance and ARPES measurements^ 
with their curves evolving smoothly across Tc, are com- 
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polling evidences in favor of a similar nature of both gap 
and pseudogap. 

Several different ideas of a superconducting phase in 
cuprates attainable by some form of percolation between 
pairs, cluster or stripes are not new^^~^^. More recently, 
a percolating approach was suggested and supported by 
the microscopic model based on a spatial dependent crit- 
ical temperature T'c(r) due to intrinsic inhomogeneities 
like pair breakers or spatial dependence of the supercon- 
ducting coupling constant proposed by Ovchinnikov et 
aP^. They have derived the density of states due to 
the spatial distribution of coupling constants. A dif- 
ferent method based on the bond percolation between 
pre-formed mesoscopic Jahn- Teller pairs was proposed 
recently^^ as a mechanism to attain a superconducting 
phase coherence . On the other hand, the spatial varia- 
tions on the LDOS due to charge disorder were taken into 
account in calculations using a generalized f — J model^^. 
In fact, new trends an ideas have been used to describe 
the effects of the intrinsic disorder and the possibility 
of a percolating phase not only for HTSC but also for 
manganites^^. 

The ideas described in the above paragraph are on 
the same lines with the scenario that we want to dis- 
cuss to explain why the occurrence of pairing and the 
phase coherent superconducting state in the high- Tc su- 
perconductors appears to occur independently: we show 
that the percolation theory, i.e., that the superconduct- 
ing state is reached through the percolation of regions 
rich in pre-formcd pairs which provides good quantitative 
agreement with the measured T* and Tc phase diagrams. 
Notice that there is a general agreement with respect to 
the experimental Tc((p)) curves, since most of them are 
obtained through the same method, namely, resistivity 
measurements. However the measured values of T* found 
in the literature are obtained through different methods 
and seems to vary considerably, depending on the spe- 
cific experimental probe used^. Such difference may be 
due to the anisotropic d-wave nature of the gap ampli- 
tude and the fact that a given experimental technique is 
sensitive to excitations at a particular wavevector mag- 
nitude. Therefore, it is natural that angle-resolved pho- 
toemission (ARPES), tunneling spectroscopy, transport 
properties such as dc resistivity and optical conductivity, 
NMR, Knight shift relaxation rate, electronic Raman, 
magnetic neutron scattering, specific heat^'^^ and recent 
vortex-like Nernst signal measurements^* yield different 
values for the onset of vanishing gap temperature T* . 
Since the origin of the charges inhomogeneities are not 
known, T* could be also dependent on the way a given 
sample is made. 

We must also emphasize that the percolating approach 
is not only suitable to yield good quantitative agreement 
with the HTSC phase diagrams but, and perhaps more 
important, it provides interesting new physical insight 
on a number of phenomena detected in these materials: 
the variation of the measured pseudogap magnitude with 
the temperature^, the decreasing of the zero temperature 



superconducting gap Aq while Tc increases for under- 
doped compounds'''^^, the downturn of the linear depen- 
dence of the resistivity with the temperature for under- 
doped samples and the increase of Hall carriers with the 
temperature, mostly measured for the optimally doped 
and underdoped compounds^°, the downturn of the lin- 
ear specific heat coefficient,.... These properties will be 
discussed in more detail in section IV. 

This paper is divid(;d as follows: in section II, wc intro- 
duce the charge distributions appropriate to mimic the 
real charge distribution inside the material and reproduce 
the T* and Tc phase diagram. In section III, we derive the 
phase diagram for Bi2212, LSCO and YBCO and make 
comparison with the experimental data. We use a mean 
field BCS like method with an extended Hubbard Hamil- 
tonian to derive the onset of vanishing gap at T*, but 
it should be emphasized that the percolating approach 
introduced here is independent of any method of calcu- 
lating the superconducting pair formation. In section IV, 
we comment on the applications and implications of the 
percolation theory to several physical properties of the 
HTSC and we finish with the conclusions in section V. 



II. THE CHARGE DISTRIBUTION 

The consciquence of the microscopic charge inhomo- 
geneities distribution in the Cu02 planes, either in a 
striped or in a granular configuration, is the existence 
of domain walls between the two phases which arc spon- 
taneously created in the planes^^^^"^: regions which are 
heavily doped or hole-rich may form the stripes and 
others regions which are hole-poor are created between 
these charge-rich stripes. The exactly form of these 
charge distributions is not well known but we can get 
some insights from a number of recent experiments. 
The neutron powder diffraction^'^'^° suggests that the 
charge inhomogeneities modifies the Cu — O bond length, 
leading to a distribution of bond lengths for optimal 
and underdoped compounds. Scanning tunneling mi- 
croscopy/spectroscopy (STM/S)^'** on optimally doped 
Bi2Sr2CaC20s+x measures nanoscale spatial variations 
in the local density of states and the superconducting 
gap at a very short length scale of « I4A. These results 
suggest that instead of a single value, the zero temper- 
ature superconducting gap assumes different values at 
different spatial locations in the crystal and their statis- 
tics yield a Gaussian distribution*. New high resolution 
STM measurements^ have revealed an interesting map 
of the superconducting gap spatial variation for under- 
doped Bi2212. Their data is compatible, as they pointed 
out, with a granular superconducting grains separated by 
non-superconducting regions. Based on their measured 
gap histogram (see Fig. 4b of Lang et al^) and on similar 
data of Pan et al*, wc can draw some insights on the real 
charge distribution inside a HTSC. 

Thus, in order to model the above experimental ob- 
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servations and to be capable of performing calculations 
which may reproduce the measured phase diagrams, 
namely T* and Tc for a given family of compounds, we 
used a combination of a Poisson and a Gaussian distribu- 
tion for the charge distribution p{r). In fact, each type 
of distribution has a convenient property which we used 
below: The width of a Gaussian distribution is easy to be 
controlled and expresses the degree of disorder but since 
it is symmetric, it would imply in a rather small width 
to the low density compounds. On the other hand, a 
Poisson distribution starts sharply and has a long tail 
which is convenient to deal with these these low densi- 
ties compounds with their experimental measured large 
degree of inhomogeneities. Thus, for a given compound 
with an average charge density {p), the hole distribution 
P(/9; (p)) or simply P{p) is a histogram of the probability 
of the local hole density p inside the sample, separated 
in two branches or domains. The low density branch 
represents the hole-poor or non-conducting regions and 
the high density one represents the hole-rich or metal- 
lic regions. As concerns the superconductivity, only the 
properties of the hole rich branch are important since the 
current flows only through the metallic region. 

Such normalized charge probability distribution may 
be given by: 

P{p) = {Pc - P) oxp[-(p - pe)V2(tT-)']/[(a_)2(2 - exp 

{-{p,f /2{a.f))] for 0<p<pe (1) 
P{p) = for pc< p< pm (2) 
P{p) = (P - Pm) exp[-(p - pra)y2{a+f]/[ia+f{2 - exp 
(-(pe)V2(fT-)'))] for p™<p (3) 

The values of cr_ (cr+)controls the width of the low 
(high) density branch. Here, pc is the end local density 
of the hole-poor branch. is the starting local density 
of the hole-rich or metallic branch. Both p^ and pm are 
shown in Fig. 2 for a particular case. Each compound will 
have a its value of p„i. Thus, a compound having regions 
with different hole concentrations which are distributed 
according to the above probability density function P{p), 
has an average value which depends on the parameters 
of the distribution. We can show that the average value 
of the charge density is 

(p) = [a+sqrt{'7r/2) + Pm + Pc- (J-sqrt{'K /2) x 
er function (pc / sqrt{-K / 2))] / 
[(a+)2(2-exp(-(pe)V2(a-)'))]. (4) 

here erfunction is the error func;tion. For most com- 
pounds, due to the small values of a+, (p) « p^, see 
Table I below. Indeed, for a compound with average 
density (p), the values of c7±are chosen in order that per- 
colation in the hole-rich branch occurs at a given density 
called Pp. 

As the doping level (p) increases, the compounds be- 
come more homogeneous and therefore, we take larger 
values of (t_ and smaller values of the metallic branch 



width (T+. We take pc ~ 0.05 as the end of the low 
branch because it corresponds to the onset of supercon- 
ductivity to mostly cuprates. However, pc ~ 0.1 for the 
Bi2212 family. During the preparation of this work, we 
discovered that the charge probability distribution for the 
metallic branch, introduced as an ansatz in the beginning 
of our calculations, but based on the results of Pan et 
al^, was confirmed by the very fine STM measurements 
of Lang et al^. They measured the gap probability dis- 
tribution which has a Gaussian form but with a long tail 
which is more characteristic of a Poisson type distribu- 
tion. Since there is not any experimental measurements 
for the hole-poor branch, fT_ is just a free parameter. 
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0.080* 


0.033 


0.050 


0.09 


0.245 




0.100 


0.034 


0.050 


0.10 


0.240 




0.120 


0.037 


0.050 


0.12 


0.238 




0.160* 


0.057 


0.040 


0.16 


0.230 




0.200 


0.078 


0.026 


0.20 


0.239 




0.220 


0.090 


0.024 


0.22 


0.245 




0.260* 


0.04(G) 


0.04(G) 


0.26 


0.268 



TABLE I. Some selected parameters for different doping 

level compounds. Notice that they are characterized by the 
average hole density (p) given in the fourth column. Notice 
also that the values of the metallic branch width a+ decreases 
with (p), showing that the degree of disorder decrease with 
the doping. The * marks the distributions plotted in Fig.2. 
The value of f7+ = 0.04 for (p) = 0.16 is taken from the 
experimental data of Pan et al®. 



4 



We can get a reliable estimation of the (t+ val- 
ues from the experimental STM/S Gaussian histogram 
distribution for the local gap^ in an optimally doped 
Bi2Sr2CaCu20s+x {Pm « 0.32). From Fig.2d of Ref.7, 
one sees that the half-width for the metallic branch is 
about half of the mean doping value, that is 2(T+ = 0.16, 
since the local gap is proportional to the local hole den- 
sity. As mentioned, the metallic branch distribution has 
exactly the same form of the STM histogram data of Lang 
et al^ which express the inhomogeneous superconducting 
regions in Bi2212. For the LSCO system, assuming that 
it has the same STM/S Gaussian histogram distribution, 
the optimally doped compound ((p) ~ 0.16) must have 
2a+ = 0.08. Since the a+ values are related with the 
degree of inhomogeneities, we estimate the a+ values for 
the other compounds using the fact that it decreases with 
the hole density. Below we show some selected distribu- 
tions for an underdoped, optimally doped and overdoped 
compounds. 

Poisson-Gaussian Distribution 
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FIG. 2. Model charge distribution for the inhomogeneities 
or stripe two phcise regions. The low density insulating (an- 
tiferromagnetic) branch is near p = 0. The high density 
hole-rich region starts at the compound average density pm 
which are indicated by long arrows for some selected com- 
pounds, pp, indicated by the short arrows, is the density 
where percolation can occur. 

The density distribution is normalized to unity. For 
compounds like LSCO with average densities of p 0.25, 
that is, in the far overdoped region, we use a single Gaus- 
sian distribution which reflects the more homogeneous 
character of these compounds. For Bi2212, since the op- 
timally doped compound is about twice that of the other 
families, the single Gaussian distribution is appropriate 
for compounds with p ^ 0.5. In the next section we 
show how to estimate as function of the local den- 
sity p{r). As already discussed, HTSC compounds are 
intrinsic inhomogeneous, and each region in the material 
may have a different superconducting critical tempera- 



ture Tc{r). The maximum of these temperatures Tc(r) is 
the pseudogap temperature since it is connected with the 
opening of the largest gap, that is, the T* of a given sam- 
ple with (p), i.e. r*((p)). Eventually for temperatures 
below T*((p)) the regions with local density between p™ 
and Pp becomes superconducting and the supcrcurrent 
may percolate through the sample. Therefore, the local 
Tc{pp) is the maximum temperature at which a dissipa- 
tionless current can flow through the system and which 
is identified as Tc{{p)), the experimental superconducting 
critical temperature. 

According to the percolation theory, the site percola- 
tion threshold occurs in a square lattice when 59% of 
the sites are filled^^ . Thus, we find the density where the 
hole-rich branch percolates integrating J P{p)dp from p^ 
till the integral reaches the value of 0.59, where we de- 
fine Pp. Below Tc{{p)) the system percolates and, conse- 
quently, it is able to hold a dissipationless supercurrent. 
We should point out that 0.59 is an appropriate value 
for a single layer cuprate. The site percolation threshold 
is 0.16 for a simple cubic lattice^^. Therefore, for a two 
layers system like Bi2212, the site percolation threshold 
may be less than 0.59, but since the true value is not 
known, we will use 0.59 as a first approximation. In the 
Table I, we show some of the parameters used for cer- 
tain sample and we plot the distribution, as discussed 
above, in Fig. 2. Notice that the Table I and Fig. 2 are 
for a HTSC system with optimum doping (p) « 0.16 like 
LSCO or YBCO. For Bi2212, the optimum hole doping is 
about twice this value^^ and in this case new parameters 
must be used. Some charge distributions used for the 
calculations with the Bi2212 family, similar to the ones 
shown in Fig. 2 above, were also derived^®. 



III. THE PHASE DIAGRAM 

The experimental based phenomenological charge dis- 
tribution probability introduced in the last section en- 
able us to estimate the appearance of pre-formed pair 
at different regions in the crystal provided that we know 
the onset of superconductivity T*(p(r)) There are sev- 
eral different approaches which can be used to obtain 
T*(p(r)). We can simply use the experimental measured 
pseudogap temperature T* = T*{(p)) for several different 
compounds through the identification T* {p{r))=T* {{p)) 
which is probably the best estimation. Another possibil- 
ity is perform a theoretical calculation as we have done in 
the past^^'"^^ for the Tc of a given compound with density 
(p) and taker*(p(r)) =Te((p)). 

Strictly speaking, due to the non-uniform charge dis- 
tribution, we do not have translational symmetry and 
we should use a method which takes the disorder into 
account^^'^^'^^'^^'^^. However such theories requires the 
introduction of a phenomenological random potential to 
simulate the disorder. Since our purpose here is to 
demonstrate that phase coherence is attained by perco- 
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lation and it is independent of the pairing mechanism, we 
will take the simplest theoretical approach''^: we use a 
BCS-type mean-field approximation in fc-space on a two- 
dimensional square lattice with uniform carrier density p 
to estimate the onset temperature of the superconducting 
d-wave gap, i.e. T*{p). Then we take T*(p(r)) = T*{p) 
in order to make a map of the superconducting region. 
We have to bear in mind that this is an approximation 
but which is worthwhile since it capable to reproduce the 
experimental T* values and is reasonable if the size of the 
grain or stripes domains are larger than the typical pair 
coherent length. In fact, the latest STM/S results yields 
grain boundaries of the order of 100-10 nanometers while 
the typical coherent lengths are of the order of angstroms. 
Thus, this approach, with appropriate choice of experi- 
mental or calculated parameters, was used before to de- 
rive the Tc{p) curves^^"^^ for some HTSC families, since 
Tc was, as in normal superconductors, taken as the onset 
of vanishing gap and p was taken as (p). A two dimen- 
sion extended Hubbard Hamiltonian in a square lattice 
has been used to model the quasi-bidimensionality of the 
carriers motion through the CUO2 planes^^~^^'*° and is 
given by 

H = - ^ tijc\„Cj^ + ^ X] '^'T'^n 

+ X ^yC-^4<^'CjCT'CiCT, (5) 

where tij is the nearest-neighbor and next-nearest- 
neighbor hopping integral between sites i and j; J7 is the 
Coulomb on-site correlated repulsion and Vij is the at- 
tractive interaction between nearest-neighbor sites i and 
j. a is the lattice parameter. 

Using the well known BCS-type mean-field approxi- 
mation to develop Eq.(5) in the momentum space, one 
obtains the self-consistent gap equation, at finite tem- 
peratures^^ 

with 

E^ = ^ei + ^i, (7) 

which contains the dispersion relation £k, and the inter- 
action potential Vkk' which comes from the transforma- 
tion to the momentum space of Eq.(5)'^^''*". In the calcu- 
lations we have used a dispersion relation derived from 
the ARPES data^^ with five neighbors hopping integrals. 
The hopping integrals could also be estimated from band 
structure calculations^^. The interaction potential may 
be given by^^.^o. 

Vkk' = U + 2V cos{kxa) cos(fc^a) -|- 2V cos{kya) cos(A;^a). 

(8) 



The substitution of Eq.(8) into Eq.(6) leads to appear- 
ance of a gap with two distinct symmetries'^^: 

Ak(T) = A(T)[cos(fca,a) ± cos(fcj,a)], (9) 

where the plus sign is for extended-s wave and the minus 
sign, for d wave symmetry. In accordance with Ref^® one 
observes that the d wave part of the gap do not depend 
on the coupling constant U , depending only on V . Here 
we deal only with the d-wave which is the more accepted 
pseudogap symmetry^. 

Using the same BCS-type mean-field approximation, 
one obtains the hole-content equation^^ 

where < p < 1. 

Eq.(lO) together with the gap equation (6) must be 
solved together self-consistcntly. They may be used to 
derive the onset of vanishing gap temperature as func- 
tion of a given value of the density of carriers p. This 
procedure was used in the past to derive, with appropri- 
ate set of parameters, the Tc{p) phase diagram for dif- 
ferent HTSC systems for a single type^^'^'''^°'^'' and for 
mixture of different order parameter symmetries'^^ . 

In the present work, we assume a different view, p here 
is not the compound average hole density but the lo- 
cal hole density p{r). For a given value of p{r), we use 
Eq.(lO) and Eq.(6) to calculate the local onset temper- 
ature of vanishing gap. It is more appropriate to call 
this temperature as T*{p); therefore, T*{p) = Tc{p) and 
hereafter we will deal only with T*{p). On the other 
hand, the pseudogap temperature of the compound is 
T*((p)) =max{rc(p)}, or simply T* . 

Below, in Fig. 3, we present the results for the 
La2-xSrxCuOi family. We have used a dispersion rela- 
tion derived from the Schabel et al^^. In their notation, 
the hopping parameters arc: t = ti=0.35eV, ^2/^1=0. 55, 
i3/ti=0.29, t4/ii=0.19, t^/ti=Qm. The magnitude of 
the attractive potential was set V/t=-QAQ in order to 
give a reasonable agreement with many measured val- 
ues of T*. As we have already mentioned, different ex- 
periments yield completely different results for T*. The 
phase diagram depend on the hopping parameters and 
on the attractive potential, which is the free parameter 
on this type of calculation. Thus, varying V/t makes the 
values of the calculated T* {p) change and the position of 
the optimal density depends on the hopping parameters. 
The hopping values used below are within the variation 
estimated by different methods (band structure, ARPES 
measurement, etc.)"*^'*^. 

The theoretical curves in Fig. 3 arc derived in the fol- 
low way: For each local density p > pm, that is, in- 
side the metallic branch, T*{p), calculated with the mean 
field equations, is a decreasing function. Thus the max- 
imum value of T*{p) is equal T*{pm). This is the onset 
temperature of the superconducting gap (in the metallic 
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branch) and therefore T*{pm) = T*{{p)). To show this, 
we draw an arrow in Fig. 3 showing how T*[{p) = 0.16) is 
obtained. The T*{{p)) curve calculated with the above 
parameter is in good qualitative agreement with the dis- 
played experimental data^^'^^. The value of supercon- 
ducting critical temperature Tc{{p)) is estimated in the 
following manner: we calculate the maximum temper- 
ature at which the superconducting region percolates 
in the metallic branch. This percolation occurs when 
all the clusters with local density between pm and pp 
are superconducting. Fig. 3a shows how this happens 
for (p) = 0.16 which yields pp = 0.23. This value of 
Pp can be seen in the panel following the arrow which 
shows that T*(0.23) is equal the superconducting criti- 
cal temperature of the compound, rc((0.16)). In general 
Tc{{p)) = T*{pp) and the system will be superconduct- 
ing when submitted to any temperature below T*(pp) 
(sec Fig. 3b). Below we show the results for the LSCO 
family. A similar curve was also studied which the T* 
is in reasonable agreement with the T^^^ extracted from 
the high-resolution dilatometry data^^ for Y Ba2Cu^Ox- 




0.05 0.1 0.15 0.2 0.25 0.3 
P.<P> 



FIG. 3. Phase diagram for the LSCO family. To explain 
how T*((p)) and Tc((p)) are obtained, we plot in (a) the prob- 
ability distribution for the optimal compound with (p) = 0.16, 
P(p, 0.16). The arrows shows T*(0.16) and the percolation 
threshold at pp = 0.23 with Tc((p> = 0.16) = Tc(0.23). The 
experimental data are taken from Ref.*** and Tonset is taken 
from the flux flow experiment of Ref.'^* (open squares). No- 
tice that values of r*((p)) are the same of Tc{p) but T* refers 
to the compound density (p) and Tc refers to the local density 
P- 



We have also used this procedure to calculate the phase 
diagram for the Bi2212 family which is in agreement with 
the experimental data as one can sec in Fig. 4. The de- 
tails of the calculations are similar to the ones described 
above and can be found in ref.19. 
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FIG. 4. a) The probability distribution for the optimal 
compound with (p) = 0.32. b)The calculated local onset of 
vanishing gap T*(p(r)) and the compound pseudogap tem- 
perature T*((p)). The thick arrow shows T*((p)) and the 
dot-dash arrow shows how Tc{{p)) is determined from T*(pp). 
The experimental points and the symbols are taking from 
Ref48. Notice that values of T*((p)) refer to the compound 
average density (p) and T* (p) refers to the local onset of van- 
ishing gap at a cluster of density p. 

It is worthwhile to mention that the percolating ap- 
proach is independent of the above mean field calcula- 
tions for T*{p). It can be used in connection with any 

method which yields a T*{p) curve. The only require- 
ment is that it would not cross the Tc{{p)) curve as some 
have proposed^^. Thus the percolating approach could 
be used with some experimental measurements or oth- 
ers theoretical methods to calculate T*{{p)), like, for in- 
stance, the calculations made with the Hubbard-Holstein 
Hamiltonian^" . 



IV. DISCUSSION 



There are several HTSC phenomena which are not well 

understood and can be explained with the percolating 
ideas, and with our model and calculations: 

1-ARPES measurements^'^^ have revealed the anoma- 
lous behavior of the zero temperature gap Ao((p)) which 
decreases steadily with the doping (p) although Tc in- 
creases by a factor of 2 for their undcrdopcd samples. 
This overall relation between Ao((p)) and Tc is totally 
unexpected since it is well known that normal supercon- 
ductors have a constant value for the ratio 2Ao/fc_Brc, 
being 3.75 for usual isotropic order parameter and 4.18 
for dx2-y2 wave solution^^. 
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FIG. 5. The zero temperature gap for 9 samples as mea- 
sured by Harris et aP* and our calculations. 

At low temperatures, since the superconducting phase 
percolates through different regions, each one has a given 
Ao(r). Thus, tunnehng and ARPES experiments detect 
the largest gap present in the compound. Consequently, 
Ao{{p)) must be correlated with the onset of vanishing 
gap T*{{p)) which is the largest superconducting temper- 
ature in the sample, and should not be correlated with 
Tc{{p)). As we show in Fig. 5, correlating the values plot- 
ted in Fig. 3b for T* ((/?)) with Ao((p)), we arc able to 
give a reasonable fit for the data^'^^ on Dy — BSCCO 
and explains the different energy scales pointed out by 
several others authors. 

2- Transport experiments have also been used to study 
the pseudogap^. The underdoped and optimum doped 
high- Tc oxides have a linear behavior for the resistiv- 
ity in the normal phase up to very high temperatures. 
However, at T* there is a deviation from the linear be- 
havior and the resistivity falls faster with decreasing 
temperature^'^*. This behavior can be understood by the 
increasing of superconducting cluster numbers and size, 
as the temperatures decreases below T* . Each supercon- 
ducting cluster produces a short circuit which decreases 
the resistivity below the linear behavior between T* and 
Tc- To obtain a quantitative fitting of this effect, we are 
presently working on a simulation for the resistivity of a 
linear metallic medium with short circuited regions which 
varies with the temperature. The linear behavior of the 
resistivity above T* may be also be explained by a per- 
colation procedure. The conductivity of inhomogeneous 
systems has been considered from different percolative 
approaches in order to obtain the correct temperature 
dependence^^'^^. The temperature functional form de- 
pends on how the percolation is achieved in the system. 

3- Several measurements made in the presence of a 

magnetic field seem to agree with the percolating sce- 
nario; The magnetotransport measurements^'^ in a 
La2-xSrxCu04 film just above the irreversibility line is 



in agreement with the notion that the magnetic field pen- 
etrates partially in the superconducting regions and de- 
stroys some of the superconducting clusters in the film. 
We have already mentioned the recent measurements of 
magnetic domains above Tc which has been interpreted as 
a diamagnetic precursor to the Meissner state, produced 
by performed pairs in underdoped La2-xSrxCuOi thin 
films^'^. The existence of superconducting clusters be- 
tween T* and Tc easily explains the appearance of local 
diamagnetic or Meissner domains and, if there is a tem- 
perature gradient in the sample, the local flux flows and 
produces the dynamic flux flow state^*. 

4- Another important consequence which follows is that 
the superconducting pairing mechanism should be more 
easily investigated by experiments performed mainly at 
T*. A such experiment was accomplished by Rubio 
Temprano et al^**, which measured a large isotope ef- 
fect associated with T* and an almost negligible iso- 
topic effect associated with T^. in the slightly underdoped 
HoBa2CiiiOs compound. The results strongly support 
the fact that electron-phonon induced effects are present 
in the superconducting mechanism associated with T*. 
Bussmann-Holdcr et al have also calculated T* as func- 
tion of a phonon induced gap^^. 

In order to gain further insight on the nature of 
the pair potential, we have measured the resistiv- 
ity under hydrostatic pressure on an optimally doped 
Hgo.82Reo.isBa2Ca2Cu30s+s^^ sample. The data in- 
dicate a linear increase of T* with the pressure at the 
same rate as Tc- In the context of our theory, this result 
might be in agreement with the phonon induced mech- 
anism: the inhomogeneities local charge densities in a 
given compound yield varying values for the Fermi level, 
broadening N{Ep)'^^ . The applied pressure on a cuprate 
with an inhomogeneous charge distribution is also ex- 
pected to broaden the density of states^^ N{Ep), and 
the main effect of an applied pressure on Tc is an in- 
crease of the phonon or Dcbye frequency. This is seen 
through the linear increase of T*^^' which also provides 
a very interesting physical explanation on the origin of 
the linear pressure induced intrinsic effect^^'^'^'^^, usually 
postulated to explain the raise of Tc above its maximum 
value. 

5- The anomalous behavior found mostly in under- 
doped compounds is related with the temperature be- 
havior of the Hall coefficient''" Rh and the Hall density 
uh which is proportional to the average hole density (p) . 
While in normal metals nn is independent of the temper- 
ature, in the normal phase of many HTSC, uh increases 
monotonically as the temperature increases and saturates 
at high temperatures^°'^^. This anomalous behavior can 
be explained by the existence of superconducting islands 
above Tc which gradually turn into metallic phase as the 
temperature is raised. Such superconducting regions can 
be regarded as empty spaces in a disordered hole-rich 
(metallic) and hole-poor (insulator) backgroimd. There- 
fore, as the metallic region increases, it increases also the 
number of carrier. However, we have to point out that 
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the saturation seems to occur near the room tempera- 
ture, which is much higher than some of the measured 
values for T* . 

6- The pseudogap has also been detected through the 
specific heat coefficient which undergoes changes below 
T*. It is well known that the electronic specific heat 
7 term of the normal phase is a material dependent 
and temperature independent constant. However, sev- 
eral measurements have detected^'^®'^^ the suppression in 
the 7 term in mostly underdoped compounds of different 
families. In the percolating approach, with the existence 
of superconducting clusters below T*, we can qualita- 
tively explain such depression. While the metallic region 
specific heat behaves as 7T, the superconducting region 
contributes with terms proportional to exp(— A(r)/T), 
which decreases much faster with T. As the temperature 
falls below T* and the superconducting region increases, 
its contribution to the specific heat becomes more impor- 
tant and produces the overall downturn in the 7 coeffi- 
cient. For overdoped samples, since T* is very close Tc, 
this effect may be difficult to be detected and this is the 
suppression of the 7 term was not seem for overdoped 
compounds^^'^^. 

7- Magnetization measurements above in oriented 
powder of Yi-xCaxBa2CusOy has revealed a pattern^° 
which cannot be explained by the usual fluctuation of 
the magnetization order parameter^^. Such pattern were 
interpreted by the formation of superconducting regions 
above the percolating threshold^" as discussed in detail 
by the percolating model above. We are presenting using 
the Ginsburg-Landau formalism®^ in connecting of the 
hole charge distribution introduced in Section II in or- 
der to provide a quantitavive calculation for such anoma- 
lous induced diamagnetism. The results will be published 
elsewhere together with a new analysis on the behavior 
of Hc2 for high- Tc superconductors. 



V. CONCLUSIONS 

We have demonstrated that the percolating approach 
for local pre-formed pairs due to an inhomogeneous 
charge distribution on the CUO2 planes provides a novel 
and general approach to the phase diagram of the high- 
Tc cuprates superconductors. The pseudogap is regarded 
as the largest superconducting gap among the supercon- 
ducting regions in an inhomogeneous compound. The 
critical temperature Tc is the maximum temperature at 
which these superconducting regions percolate. We have 
shown, through the calculations presented in this paper, 
that this approach is suitable to reproduce the measured 
Tc and T* phase diagrams for several cuprates using a 
phenomenological real charge distribution drawn from 
the STM/S data. 

One of the advantages of our approach is that the pair 
formation and the phase coherent superconducting state 
are studied independently. The phase coherence is at- 



tained by percolation. On the other hand, the pair for- 
mation was studied by a mean field BCS like method in 
order to estimate the onset of superconducting gap T* 
for any doping level. This was the simplest method and 
the calculations could have been done with other theoret- 
ical approaches. However, notice that the pair formation 
mechanism is independent of how the coherent percola- 
tive phase is reached. This general approach is in agree- 
ment with new trends and ideas^*''^^, which are based on 
new data and new experimental facts that mostly HTSC 
are inhomogeneous materials. 

The method developed provides also new insights and 
introduces a new interpretation on several typical anoma- 
lous properties of HTSC: like the dependence of the zero 
temperature gap Ao(r) with the hole concentration, the 
downturn of the linear dependence to the resistivity with 
the temperature, the linear dependence of the hole con- 
centration with the temperature, the suppression of the 
specific heat coefficient 7, the dependence of T* with the 
pressure, and others new implications which are presently 
being studied and which will be discussed in the future. 
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